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£S| ' Abstract 

a' 



This paper focuses on the analysis of average Gaussian error probabilities in certain fading channels, 
i.e. we are interested in E[Q(y/Jry] where Q(-) is the Gaussian Q-function, p is a positive real number 
and 7 is a nonnegative random variable. We present a unified analysis of the average Gaussian error 
probability, derive a compact expression in terms of the Lauricella function that is applicable to a 



broad class of fading channels, and discuss the relation of this expression and expressions of this type 
recently appeared in literature. As an intermediate step in our derivations, we also obtain a compact 



expression for the outage probability of the same class of fading channels. Finally, we show how this 
unified analysis allows us to obtain novel performance analytical results. 
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■ Index Terms 

Performance analysis, average bit error probability, outage probability, fading channels, Lauricella 
functions. 
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i_l ', I. Introduction 

a 

For many decades, communication theorists have analyzed the performance of single-channel 
and multichannel receivers in a fading environment. Among performance measures, the average 
bit error probability (BEP) is perhaps the one that is most revealing about the nature of the 
system behavior. The average BEP is defined as the average of the conditional BEP over the 
fading statistics. The conditional BEP of a fading channel is often equivalent to the BEP of an 
additive white Gaussian noise (AWGN) channel, e.g. systems employing M-AM or QAM under 
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ideal coherent detection [1]. In such cases, the average BEP can be computed from the statistical 
expectation of Gaussian error probabilities with respect to the fading distribution. 

This paper focuses on the analysis of average Gaussian error probabilities in fading channels, 
i.e. we are interested in E[Q(y/py] where Q(-) is the Gaussian Q-function, p is a positive 
real number and 7 is a nonnegative random variable. We may note that the final average BEP 
expression is usually expressed as a weighted sum of average Gaussian error probabilities, e.g. see 
[2] for QAM with Gray mapping; however, for clarity these details will be overlooked here. The 
literature concerning average Gaussian error probability calculations is now quite voluminous. 
The most complete account of this problem is found in [1]; nevertheless, a great deal of new 
results have appeared after the publication of [1]. Some of these new closed-form results involve 
the Lauricella function \ e.g. [3]-[5]. The approach identified in these new contributions is 
based on applying an appropriate change of variable to the expression obtained by the well- 
known moment generating function (MGF) method developed in [1]. Therefore, these results 
raise the natural question of whether the involved fading distributions share a common property 
that leads to this particular mathematical form. 

In this paper we identify the common property of certain fading distributions that leads to 

in) 

average Gaussian probabilities expressed by the Lauricella function F D '. Such property is related 
to the form of the associated MGF. We derive unified expressions for the average Gaussian error 
probability and the outage probability, which are applicable to a large class of fading distributions. 
Moreover, this unified analysis provides a systematic method for obtaining new analytical results. 

The remainder of this paper is organized as follows. The unified analysis is presented in 
Section II. In Section III we apply our analysis to derive both published and novel analytical 
results. Finally, some conclusions are given in Section IV. 

II. Unified Analysis 

In this section we derive the key results of this work. Some comments on notation are in 
order. For an arbitrary function <\>[x) we denote the Laplace transform as C[(f)(x); s]. As in [1], 
we define the MGF of a nonnegative random variable 7 as M 7 (s) = E[e S7 ] = £[/ 7 (7); — s], 
where s 6 C and / 7 (7) is the probability density function (PDF) of 7. 

The following definition will be very useful for our purposes 1 . 

'This terminology is inspired by the geometric programming theory. 
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Definition 1 (Monomial and Posynomial MGF): A nonnegative random variable 7 has a posynomial 
MGF if its MGF has the form 

K n k / \ -b k ,i 

M,(s) = J2^H 1 , (1) 

k=l i=l ' ' 

where the involved parameters satisfy the following conditions 

Re[a kti ] > 0, % = 1, 2, . . . , n k and k = 1, 2, . . . , K, 

£>,i>0, fc = l,2,...,X, 
' »=i ^ 
if 

, fc=i 

For the special case K — 1, the MGF will be called monomial. 

In the interest of brevity, we will say a random variable or a distribution is monomial or 
posynomial if its associated MGF is monomial or posynomial, respectively. In addition, the 
coefficients in Definition 1 will be referred as characteristic coefficients and the conditions given 
in (2) will be called compatibility conditions. The necessity of these conditions will be revealed 
along this section 2 . As shown in Table 1, identifying 7 with the instantaneous signal-to-noise 
ratio (SNR), we observe that many common fading distribution have monomial MGF. 

To grasp the generality of the results derived in this section, we may define the average 
Gaussian error probability in a more fundamental form. Integral transforms theory provides 
a convenient framework for our purposes. In [11] a new integral transform, called Shannon 
transform was considered in relation to the ergodic capacity analysis. Regarding the problem 
treated here, we may consider the following integral transform. 

Definition 2 (Gaussian Q-transform): Given a nonnegative absolutely continuous random vari- 
able 7 with PDF / 7 (7), we define 

Q>) = Q[/ 7 ( 7 );p]=e[Q(Vp7)], (3) 

where p is a real nonnegative number. 

Some straightforward properties of the Gaussian Q-transform, including an inversion formula, 
are collected in the following result. 

2 We anticipate that the first and second compatibility conditions will be used in the proof of Theorem 1, while the third 
condition is necessary to satisfy M~/(0) = 1. 
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Lemma 1: Let us consider a nonnegative absolutely continuous random variable 7. Then, the 
corresponding Gaussian Q-transform is always defined in the set R + = [0, 00), and the following 
properties hold: 

(i) Q-y(p) is a continuous and decreasing function in IR + such that <2 7 (0) = \ and lim Q 1 (p) = 0. 

(ii) If Q 7 (p) = $(p) then 



for any 5 > 0. 

Proof: See Appendix I. ■ 
Although the Gaussian Q-transform is an attractive concept, it is beyond the scope of this paper 
to develop these matters further. 

Under the previous theoretical framework, one can think in Gaussian Q-transform pairs. The 
key result of this work is conceptually summarized as follows: monomial distributions and certain 
expressions involving the Lauricella function are connected by the Gaussian Q-transform. 
The extension of this idea to posynomial distributions is straightforward. From an operational 
point of view, once the underlaying fading distribution is identified as posynomial, obtaining the 
Gaussian Q-transform reduces to extracting the characteristic coefficients from the MGF. The 
mathematically precise statements are given below. 

It is most common to find monomial distributions in practical problems, thus, we start pre- 
senting the results for the monomial case. 

Theorem 1: Let 7 be a monomial random variable, i.e. 



^)=n( 1 -^) > w 



where {aj}™ =1 and {&j}" =1 satisfy the compatibility conditions given in (2). Then 
(i) The cumulative distribution function (CDF) of 7 is expressed as 



ft {a t f 



R < (7) = { r 1 n x \ [ K ■ ■ ■ , b n ; 1 + b " ~ a ^> • • • > ) > ( 5 ) 



i=i 



(n) 

where T is the gamma function and $2 is the confluent Lauricella function defined in [12]-[13]. 
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(ii) The Gaussian Q-transform of 7 is given by 

drip) 




(6) 



x 



\ i=\ i=l ^ 



2Qn 
P 



where F D is the Lauricella function defined in [12] -[13]. 

Proof: See Appendix EL 
Extending this last result to posynomial random variables is straightforward. 

Corollary 1: If 7 is a posynomial random variable such that its MGF is given by (1) then 
(i) The CDF of 7 is expressed as 



n K*r 

*=i r i + £& fc)i 



b k,: 



1=1 



(7) 



i=l 



(ii) The Gaussian Q-transform of 7 is given by 



n k 

(«i -\ 1 r ( \ + J2 b k ,i 



1 " 

V fc=1 I j =1 



2 

i=l 



r ( 1 + E ^ 

i=l 




(8) 



2 + y^^fc,i,^fc,i, • • • ; &fc,ri fe ; 1 + h,f, 
i=i i=i 



2afc 1 2afc. 



, . . . , 



P P 

Proof: Repeat the steps of the proof of Theorem 1 . ■ 
Theorem 1 and Corollary 1 have important theoretical and practical consequences. Next result 

provides an integral representation for the Gaussian Q-transform of posynomial random variables 

that is convenient for numerical evaluation. 

Corollary 2: If 7 is a posynomial random variable such that its MGF is given by (1) then 
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if* 

K ( n k -\ . > J2 b k ,i 



(^) 

du. 



Proof: Check that the compatibility conditions allow us to use the Euler-type integral 
representation given in [12, p. 283, eq. 34] for the Lauricella functions in (8). ■ 

From this last Corollary it is straightforward to obtain an asymptotic approximation for 
posynomial random variables. 

Corollary 3: Let 7 be a posynomial random variable such that its MGF is given by (1). Then 



K ( n k 



v *-i J rh + Ei>w)p 



.» I r)»=i 

V i=l 

for large values of p. 

Proof: Substitute the integrand in (9) by its asymptotic approximation and use the basic 
properties of the beta function [15, p. 898-899]. ■ 
As expected, from (10) it is inferred that the underlaying diversity order of Q 7 (p) is 

p^oo log (j)) 

assuming 7 is a posynomial random variable 



n k 



liiu = min ^ > , (11) 



i=i 



III. Applications 

The mathematical tools developed in previous section provide us a unified analytical framework 
for many existing results in literature; in addition, they systematically allow us to obtain new 
analytical results. 

To calculate the average BEP, its asymptotic approximation and the outage probability we can 
follow four steps: 

1) Check if the underlaying distribution is posynomial. In such case extract the characteristic 
coefficients. 

2) Use the Theorem and Corollaries of previous section to obtain analytical expressions for 
the Gaussian Q-Transform, its asymptotic approximation and the outage probability. 

January 13, 2013 DRAFT 



3) If possible, reduce the Lauricella functions and to simpler functions. 

4) Use known formulas to express the average BEP in terms of the Gaussian Q-transform, 
e.g. see [2] for QAM with Gray mapping. 

Next, we provide some examples to illustrate how this approach allows us to derive both 
published and novel results. For brevity we will mainly focus on the first step described above. 

Example 1: Nakagami-m fading. This a well-known example [ 1 ] [ 1 4] . From Table I, we can 
extract the characteristic parameters for the Nakagami-m distribution; specifically, ai = m/7 
and bi = m. From (5) we obtain 



in 



Then, after considering [15, eq. 8.351-2], we derive the well-known formula 

r (rn, ^7) 



^7(7) = 1- V? \ • (13) 

r (m) 

Applying (6) and considering [15, eq. 9.131-1] yields 

where c(p) = p r y/(2m). This expression is very similar to that given [1, eq. 5.A.2][14, eq. A. 8], 
which was derived by a different approach. Interestingly, it has been checked numerically that 
both are equivalent; however, the author has not been able to find the transformation between 
Gauss hypergeometric functions 2 F± connecting (14) and [1, eq. 5. A. 2]. □ 
Example 2: Maximal ratio combining (MRC) with independent but nonidentically distributed 
branches. A formula for the Gaussian Q-transform over Hoyt fading channels is derived in [4] 

in) 

in terms of the Lauricella Fp function. Now we show how to perform a unified analysis for 
MRC including Hoyt distributed branches. Let us consider receive MRC with L independent but 
nonidentical distributed branches. We assume that L\ branches exhibit Nakagami-g (Hoyt) fading 
and the remainder of the L — Li branches are better modelled by the Nakagami-m fading model. 
The Nakagami-m fading model is used for both line-of-sight (LOS) and NLOS scenarios, while 
the Nakagami-g distribution is an alternative to model NLOS channels. We assume a general 
scenario where each branch has arbitrary parameters, i.e. {qj, 7j},=i for the first L 1 branches and 
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{rrij, 7j}j=L 1+ i for the remainder of the branches. Analyzing such scenario is straightforward 
with the results derived in this work, after observing that the associated MGF is monomial 

Li , 9 2 v -1/2 2Li , v -1/2 Li+L / - x-mi 

= n (i - i^p) n (i - rh'-v n - £«) • <*> 

i=l v Hl 7 i=Li+l v 1 / j=2Li+l v * 7 

This scenario can be further generalized to consider large-scale time variations. Let us assume that 
the number of branches with Nakagami-g fading L x is also random. Then, if Pr[/] represents the 
probability of having / branches with Nakagami-g fading, the associated MGF is now posynomial 

L If .9 \ -1/2 




i=i 




(16) 



2 -%A n fi-^ 



-mi: 



A system employing MRC with independent 77-yu. distributed branches has been recently analyzed 
in [5]. This system generalizes the Nakagami-g/Nakagami-m scenario described above and it 
requires to estimate the set of 3L parameters {rjj, fij, 7j}j=i- Since the Nakagami-g/Nakagami-m 
scenario requires 2L channel parameters, it is a reasonable alternative for certain applications. 
Interestingly, if we apply Theorem 1 to the monomial MGF considered in [5, eq. 2], we can 
extend the performance analysis in [5] with a closed-form expression for the outage probability 
and an asymptotic approximation for the average BEP. □ 
Example 3: Orthogonal space-time block codes (OSTBC) in spatially correlated multiple- 
input multiple-output (MIMO) channels. In [16] the underlying MGF was derived for different 
systems using OSTBC in shadowed Rician MIMO fading channels. Such MGFs are posynomial, 
thus, all results derived in previous section are applicable, leading to novel analytical results. To 
illustrate this idea, we consider the particular MGF for correlated LOS component and spatially 
white scattered component [16, eq. 23] 

ntn r 



v ' 1=1 x ' 



ntn r 



(17) 



(1 + as y^n t -Q (1 + (fl + b)s ym JJ (1 + as) ™ 



i=l 



where n t ,n r ,a,b and m are channel and system parameters defined in [16]. After rewriting the 
MGF given in [16, eq. 23], we clearly observe in (17) that it is a monomial MGF. □ 
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IV. Conclusions 

A solid analytical framework for the computation of average Gaussian error probabilities has 
been derived in this paper. This quantity allows us to obtain the average BEP in a variety of 
wireless communications systems. It has been shown that if the MGF of the underlaying fading 
distribution has certain posynomial structure, the average Gaussian error probabilities and the 
outage probability can be expressed in terms of Lauricella functions. Finally, this analytical 
framework has been used to derive known and novel results in a simple and systematic way. 

Appendix A 
Proof of Lemma 1 

The existence of the Gaussian Q-transform is clear after taking into account that Q{^/p^j) < \ 
for every p > 0. 

(i) Continuity follows form the fact that the kernel Q(^yp ; y) is continuous and bounded in R + . 
If pi < p 2 then Q(y/pTy) > Q(^/pla) f° r ai l 7 m thus, Q T (pi) > Q 7 (p2)- The property 
lim <2 7 (p) = follows after applying the monotone convergence theorem. 

(ii) After integrating by parts Q y (p) we can write 

$(p) = f°° e-^FJt)t-^ 2 dt. (18) 

The integral in (18) is recognized as a Laplace transform converging for 5R[p] > 0. Thus, after 
applying the inverse Laplace transform to this equality, the desired result is obtained. 

Appendix B 
Proof of Theorem 1 

bi 



(i) Since C[f 7 (t) ; s] — Y[ 1 1 + f: ) > we can write 

i=i * 1 ' 




for 3ft [s] > 0. Taking into account the first and second compatibility conditions given in (2) and 
identifying [13, p. 222, eq. 5] from (19), the expression (5) is obtained. 
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(ii) Introducing (5) in (18) yields 

{n (^o^ 1 
l= ) T x r e-iH 1 - 1 '^**™ (b 1 ,...,b n ,j2bi-,-ait,...,-a n t)dt. 
r EMj io V *=i / 

(20) 

Taking into account the first and second compatibility conditions given in (2) and identifying 
[12, p. 286, eq. 43] from (20), the expression (6) is obtained. 
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TABLE I 
Some Monomial MGFs. 



Fading Distribution 



MGF M 7 (s) 



References 



Rayleigh 



(1-7*)" 1 ; 

0^7- 



[6],[1] 



Nakagami-g (Hoyt) 



2q 2 



1 + q ilS ) y 1 + q 2 

< 7,0 < q < 1. 



-1/2 



1 



J 8 



-1/2 



[7],[l] 



Nakagami-m 



(l-lsX 
V to / 



\ ^ to. 



[8],[1] 



Rician Shadowed* 



7 



m— 1 



l-i 



A' 



7 



1 + K J \ \ mj 1 + A 

7, sC A, < m. 



[9],[1] 



1 - s- 



7 



^/2 



1 - S 



7 



n(h-H) 



-n/2 



77-/Z Physical Model 



n(h + H) 
0<7,n = l,2,... 

Format 1 =► = 2 + ^ 1 + T \ H = < v < 00 



Format 2^ h 



,H 



1 — r? 2 ' 1 — ?7 : 



:,-!<»?< 1 



[10] 



* Note that A = using the same notation as in [9],[1]. 
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